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t-h ■ Abstract 

in ; 

^ ! Suzaku observation of an ultraluminous X-ray source, NGC 2403 Source 3, per- 

formed on 2006 March 16-17, is reported. The Suzaku XIS spectrum of Source 3 was 
described with a multi-color black-body-like emission from an optically thick accretion 
disk. The innermost temperature and radius of the accretion disk was measured to be 
T in = l-08±g-g| keV and R in = 122.lj£| a 1 / 2 km, respectively, where a = (cos60°/cosi) 
with i being the disk inclination. The bolometric luminosity of the source was esti- 
mated to be Lboi = 1-82 x 10 39 a ergs s _1 . Archival Chandra and XMM-Newton data 
of the source were analyzed for long-term spectral variations. In almost all observa- 
tions, the source showed multi-color black-body-like X-ray spectra with parameters 
similar to those in the Suzaku observation. In only one Chandra observation, how- 
ever, Source 3 was found to exhibit a power-law-like spectrum, with a photon index of 
T = 2.37 ±0.08, when it was fainter by about ~ 15% than in the Suzaku observation. 
The spectral behavior is naturally explained in terms of a transition between the 
slim disk state and the "very high" states, both found in Galactic black hole binaries 
when their luminosity approach the Eddington limit. These results are utilized to 



1 



argue that ultraluminous X-ray sources generally have significantly higher black-hole 
masses than ordinary stellar-mass black holes. 

1. Introduction 

Since the era of the Einstein observatory (Fabbiano & Trinchieri 1987), luminous point- 
like X-ray sources, of which X-ray luminosities exceed a few times 10 39 ergs s _1 , were frequently 
found at off-center regions of nearby normal galaxies. These X-ray sources are called ultralu- 
minous X-ray Sources (ULXs; Makishima et al. 2000). The nature of ULXs has remained one 
of the important unresolved issues in X-ray astrophysics, for nearly three decades. 

Both observational (e.g., Mizuno et al. 2001; Kubota & Makishima 2004; Tsunoda et al. 
2006; Mizuno et al. 2007; Isobe et al. 2008) and theoretical (e.g., Watarai et al. 2000) studies 
suggest that the X-ray properties of ULXs resemble those of Galactic black hole binaries (BHBs) 
in high accretion rates. However, we currently have two distinct interpretations of ULXs; 
intermediate mass BHBs with a mass of M> 10M Q (where M & is the solar mass) radiating at 
sub- or trans-Eddington luminosities, and stellar-mass BHBs shining at highly super-Eddington 
luminosities. We urgently need to distinguish the two alternative. 

X-ray spectra of well-studied ULXs are approximated by either a multi-color disk (MCD; 
Mitshuda et al. 1984; Makishima et al. 1986) model or a power-law (PL) one (or their com- 
bination). Furthermore, several ULXs are reported to show a spectral transition between the 
MCD-like and PL-like states (e.g., Kubota et al. 2001). These properties could be explained, 
in the simplest manner, by regarding the MCD-like and PL-like states as analogous to the 
classical high/soft and low/hard states of Galactic BHBs, respectively. However, this naive 
interpretation has encountered several difficulties (e.g., Makishima et al. 2000; Mizuno et al. 
2001; Kubota et al. 2002). 

Recent observational studies (e.g, Tsunoda et al. 2006) revealed that the X-ray spectra of 
the MCD-like ULX requires a temperature profile in the accretion disk that is flatter than those 
in the standard accretion disk (Shakura & Sunyaev 1973) assumed in the MCD model. The 
flat temperature profile is theoretically predicted by a concept of slim accretion disk (Watarai 
et al. 2000), in which an optically-thick advection and/or photon trapping (e.g., Ohsuga et al. 
2005) become important due to high accretion rates. The slim disk interpretation can solve 
several basic problems with the MCD model, including too-high disk temperatures and variable 
innermost disk radii (Mizuno et al. 2001; Watarai et al. 2000) . Therefore, the MCD-like ULX is 
thought to harbor a slim accretion disk rather than a standard disk. On the other hand, photon 
indices of the PL-like ULXs (r > 2.0; e.g., Kubota et al. 2002) are steeper than those of the 
classical low/hard state BHBs (r = 1.5 - 2.0; e.g, Tanaka 1997). Because such a steep photon 
index is observed from Galactic BHBs when they are in the very high state (VHS; Miyamoto 
et al. 1991) in which Comptonization plays an important role (Kubota & Makishima 2004), 
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the PL-like ULXs are thought to be in the VHS (Kubota et al. 2002). In the recent Suzaku 
(Mitsuda et al. 2007) observation of a nearby normal galaxy NGC 1313, Mizuno et al. (2007) 
have reinforced the VHS and slim disk state interpretations of two ULXs therein, XI and X2, 
respectively. 

So far, these typical ULX properties have been observed from objects with very high 
luminosities (e.g., ^ 10 ergs s" 1 ), which require either a high BH mass or an extreme super- 
Eddington condition (or both). If, however, such ULX properties are confirmed in less luminous 
objects with a luminosity within a few times the Eddington limit of Galactic stellar-mass BHBs, 
we will become sure that the ULX phenomena can appear without highly-supper Eddington 
luminosities. Particularly important is to search for spectral transitions among such objects. 

Source 3 (hereafter, Src 3) in the nearby spiral (SABcd) galaxy NGC 2403 is one of 
the most suitable targets for our purpose. Originally discovered by the Einstein observatory 
(Fabbiano & Trinchieri 1987), it was later revealed with ASCA to exhibit an MCD-type spec- 
trum with an innermost disk temperature of T in = 1.10io]o9 keV and an innermost disk radius of 
i?i n = 130 km (Kotoku et al. 2000). Furthermore, its bolometric luminosity, L^ \ ~ 2 x 10 39 ergs 
s _1 , places it in between luminous ULX and ordinary stellar-mass BHBs. In order to examine 
whether this object exhibits the ULX-like behavior, we conducted a Suzaku observation of NGC 
2403; then the object was again found in the MCD-like state. We also utilized Chandra and 
XMM-Newton archival data of Src 3, and discovered a spectral transition that resembles those 
observed in more luminous ULXs. In the present paper, these results are utilized to strengthen 
the classification of Src 3 as a ULX with a rather modest luminosity. We employ the distance 
to NGC 2403 of 3.2 ±0.4 Mpc (Freedman & Madore 1998), as determined from observations of 
Cepheid variables. 

2. Suzaku Observation and Data Reduction 

The Suzaku observation of NGC 2403 was conducted on 2006 March 16-17, during the 
science working group phase. The X-ray Imaging spectrometer (XIS; Koyama et al. 2007) and 
the Hard X-ray Detector (HXD; Takahashi et al. 2007; Kokubun et al. 2007) onboard Suzaku 
were operated in the normal clocking mode without any window option, and in the normal 
mode, respectively. The nucleus of the NGC 2403 galaxy was placed at the XIS nominal 
position (Serlemitsos et al. 2007). 

Utilizing the HEADAS 6.5.1 software package, we reprocessed the data from the Revision 
2 processing, and created new cleaned event files. The present paper concentrates on the XIS 
data, because no significant X-ray signals were detected with the HXD in this observation. 
We screened the data under the following criteria; the spacecraft is outside the south Atlantic 
anomaly (SAA), the time after an exit from the SAA is larger than 436 s, the geometric cut-off 
rigidity is higher than 6 GV, the source elevation above the rim of bright and night Earth is 
higher than 20° and 5°, respectively, and the XIS data are free from telemetry saturation. As 
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a result, we have obtained 62.7 ks of good exposure. In the scientific analysis below, we utilize 
only those events with a grade of 0, 2, 3, 4, or 6. 

3. Results 

3.1. X-ray Image 

Figure 1 shows a 0.5 - 10 keV Suzaku XIS image of NGC 2403, superposed on an optical 
image taken from the Digitized Sky Survey (DSS; Lasker et al. 1990). Two bright point-like 
X-ray sources are clearly seen in the figure; the brighter and fainter of them correspond to Src 
3 and Source 5 (hereafter Src 5; Fabbiano & Trinchieri 1987), respectively. 

In addition to the point sources, we clearly observe faint X-ray emission, extending 
along the optical galaxy. The detailed properties of this extended (possibly diffuse) X-ray 
emission will be reported in a separate paper (Senda et al. 2008). In the summer of 2004, a 
supernova, SN2004dj, exploded in this galaxy (Nakano et al. 2004). Although X-ray emission 
from SN2004dj was detected with Chandra in 2004 (Pooley & Lewin 2004) at a 0.5 - 8 keV 
flux of 1.2 x 10~ 13 ergs s~ x cm -2 , we found no significant X-ray enhancement at the location of 
the supernova (the filled star in figure 1) above a 0.5 - 8 keV flux upper limit of 1 x 10~ 14 ergs 

— 1 —9 

s cm . 

3.2. X-ray spectrum of NGC 2403 Src 3 

We integrated X-ray signals of Src 3 within a circle denoted as Src3 in figure 1. The 
background spectrum was derived within a circle BGD3. In order to reduce possible contribu- 
tions from the extended X-ray emission, we adopted a radius of 2' (corresponding to 1.86 kpc 
at the distance of 3.2 Mpc) for these regions, which is slightly smaller than a typical integration 
region used in the Suzaku analysis of point sources (~3'). 

Figure 2 shows background-subtracted two-band X-ray lightcurves of Src 3, obtained 
with the three front-illuminated (FI) CCD cameras (XIS 0, 2 and 3; Koyama et al. 2007) of the 
XIS. Summed over the three FI CCDs, the time-averaged count rate of Src 3 is measured to be 
0.112 ±0.002 cts s" 1 and 0.078 ±0.002 cts s" 1 in the soft (0.5 - 2 keV) and hard (2 - 10 keV) 
bands, respectively. The data suggest a variation on a time-scale of > 10 ks in the middle of 
the observation, it is statistically insignificant (x 2 /d.o.f. = 10.2/14 and 11.2/14 for the soft and 
hard band, respectively). As a result, the hardness is found to be rather constant throughout 
the exposure. Below, we hence analyze the time-averaged spectrum of Src 3. 

Figure 3 shows background-subtracted XIS spectra of NGC 2403 Src 3, presented without 
removing the instrumental response. The X-ray spectra obtained with FI and background- 
illuminated (BI) CCD camera (XIS 1; Koyama et al. 2007) were both binned into energy 
intervals each with at least 100 events. Thus, the X-ray signals were significantly detected over 
the 0.4-10 keV range. The spectrum appears to be rather featureless, without any clear sign 
of emission or absorption features. 
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To analyze the spectra, we calculated a response matrix function (rmf ) and an auxiliary 
response file (arf), using xisrmf gen and xissimarf gen (Ishisaki et al. 2007), respectively. The 
reduction of the effective area due to contaminants on the optical blocking filters of the XIS was 
taken into account by xissimarf gen. We performed the spectral fitting using XSPEC version 
12.4.0. 

We fitted the Src 3 spectra with the MCD and PL models, as a standard way in analyzing 
the spectra of Galactic BHBs and ULXs. The FI and BI spectra were jointly fitted, with the 
model normalization allowed to differ between them. We found that the fluxes from BI and FI 
spectra agreed with each other, within 4%. 

Conducting the spectral fitting, we found the absorption column density, which was left 
free, to fall below the Galactic line-of-sight value of A H = 4.1 x 10 20 cm" 2 toward NGC 2403 
(Dickey & Lockman 1990). This effect may be ascribed to residual soft diffuse emission within 
the Src 3 region, which remains even after subtracting signals from the BGD3 region, due to 
brightness gradients. According to Senda et al. (2008), the diffuse emission in NGC 2403 can be 
described by thermal plasma emission model with two temperatures of ~ 0.26 keV and ~ 0.74 
keV, having a metalicity of ~ 0.44 solar. We hence took into account the contamination, by 
adding two APEC components in XSPEC with the temperature fixed at kT = 0.26 keV and 
0.74 keV, both of which were subjected to the Galactic absorption. Their relative normalization 
was fixed at the value determined with the total diffuse emission (Senda et al. 2008), while their 
summed luminosity was left free. We left free the absorption to the ULX component. 

The spectral parameters obtained in this analysis are summarised in table 1. The PL 
model was rejected with x 2 — 349.4/174, because the observed X-ray spectrum has a more con- 
vex shape, as revealed by residuals in figure 3. On the other hand, the MCD model successfully 
described the X-ray spectrum of Src 3 with \ 2 = 182.0/174. The derived absorption column 
density for the ULX component, Ah = 1. 401q!42 x 10 21 cm -2 , is higher by a factor of about 3 
than the Galactic value, but this amount of excess absorption is attributed reasonably to that 
within the NGC 2403 galaxy, and/or an additional absorber localized around the source. The 
thermal plasma components were found to have only a negligible contribution to the observed 
spectrum, with a 0.5 - 10 keV luminosity of L t h = 3.5 x 10 37 ergs s _1 . This luminosity is 
reasonable in view of the results by Senda et al. (2008). 

The innermost disk temperature of Src 3 was obtained as T m = 1.08io!o3 keV. The X-ray 
flux was measured to be 1.1 x 10 -12 ergs cm -2 s -1 in the 0.7 - 7 keV range, without correcting 
for absorption. These give an absorption-corrected bolometric luminosity of Lboi = 1-82 x 10 39 a 
ergs s _1 , with a = (cos 60°/ cosi) where % is the inclination of the accretion disk to our line of 
sight. The value of T in agrees with that obtained with ASCA in 1997 (Kotoku et al. 2000) within 
statistical errors, while Lboi is lower by about 20%. Through a relation Lboi = 47r<7r 2 n T^, with 
a being the Stefan-Boltzmann constant, the apparent innermost disk radius is calculated as 
ri n = 102.615^ km. The true innermost radius is determined as R m = 122.ll 6 ; 8 a 1//2 km, applying 
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the correction as R m = C,K, 2 r- m , where £ = 0.412 is a correction factor for the inner boundary 
condition (Kubota et al. 1998), and k = 1.7 is a spectral hardening factor (Shimura & Takahara 
1995) which represents the ratio of the color temperature to the effective temperature. 

We also examined the "variable-p" disk model, in which the dependence of the local 
temperature on the radius r from the BH is assumed to scale as r~ p , with the index p being 
a positive free parameter (Mineshige et al. 1994). While this model with p = 0.75 reduces to 
a simple MCD model, that with smaller value of p, down to 0.5, is considered to approximate 
the X-ray spectra from a slim disk (Watarai et al. 2000; Kubota & Makishima 2004; Mizuno 
et al. 2001; Mizuno et al. 2007). The fit was acceptable, but with no significant improvement 
(x 2 = 181.9/173) over the MCD fit. All the parameters became consistent with those of the 
MCD model; p = 0.73±g;^, T in = 1.09 ± 0.07 keV, and R in = 116.4±i|;| km. 

4. Analysis of Chandra and XMM-Newton Data 

In order to better understand the nature of Src 3, its intensity-correlated spectral vari- 
ations are considered to be of particular importance (e.g, Mizuno et al. 2001). Therefore, we 
analyzed the archival Chandra and XMM-Newton data of NGC 2403. 

4-1- Chandra observations 

So far, 5 Chandra ACIS exposures toward NGC 2403 have been conducted, 4 of which 
were motivated by the discovery of SN2004dj (Nakano et al. 2004). However, NGC 2403 Src 3 
was within the ACIS field of view on only three occasions. Furthermore, in the first observation 
performed on 2001 April 17 (ObsID = 2014), the ACIS spectrum of NGC 2403 Src 3 was severely 
distorted by event pileup (> 30%; Schlegel & Pannuti 2003). Therefore, we analyzed only the 
data in the remaining two observations (ObsID = 4628 and 4630), wherein Src 3 was placed 
at a relatively large distance from the ACIS aimpoint (see A6 in table 2), so that a broadened 
point spread function made the event pileup insignificant. The log of these two observations is 
shown in table 2. The data were read out in the timed-exposure mode with the standard frame 
time of 3.2 s, and telemetered in the FAINT format. In ObsID = 4628 and 4630, Src 3 was 
placed on ACIS-S3 (a BI chip) and S2 (an FI one), respectively. 

Utilizing the CIAO 3.4 software and referring to CALDB 3.3.0.1, we repro- 
cessed all the data to create new level-2 event files in the standard manner. We re- 
moved acis_detect_af terglow correction, and applied a "new" bad pixel file created by 
acis_run_hotpix l . After removing X-ray point sources detected by wavdetect within the 
ACIS field of view, we produced a 0.3 - 12 keV lightcurve covering the entire CCD chip, and 
then discarded those time intervals where the count rate exceeds 120% of that averaged during 
the observation. This procedure gave us the good exposures as listed in table 2. The ACIS 
spectra of NGC 2403 Src 3 were extracted within a circle centered on the X-ray peak, and 

1 http:/ /asc. harvard. edu/ciao3.4/guidcs/acis_data. html 
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those of BGD from a concentric annulus. The radii of these regions are also shown in table 
2. The rmf and arf files were generated, assuming an X-ray point source, with the CIAO tools 
mkacisrmf and mkarf , respectively. 

4-2. XMM-Newton observations 

So far, XMM-Newton observations of NGC 2403 have been performed 5 times. However, 
only the data obtained in three of them are public at present, of which the observational log is 
given in table 3. Although the results in the individual observations were separately reported by 
several authors (Feng & Kaaret 2005; Stobbart et al. 2006; Winter et al. 2006), we re-analyzed 
all the data in a sytematic way. The third observation (ObsID = 0164560901) was motivated 
by the explosion of SN2004dj. On all these occasions, the EPIC MOS and pn cameras were 
operated in the nominal full frame mode. The medium optical-blocking filter was adopted in 
the first and second observations (ObsID = 0150651101 and 0150651201, respectively), while 
the thick one was utilized in the third observation. In the first and second observations, we do 
not use the pn data, because Src 3 fell on the gap of the pn chips. 

We used the Science Analysis System (SAS) version 7.0.0 to reduce the XMM-Newton 
data. All data were reprocessed by emchain or epchain, based on the Current Calibration 
Files (CCF), which was latest at the beginning of our analysis (2007 July). After Read & 
Ponman (2003), we rejected high-background periods, using the point-source removed lightcurve 
with PATTERN == and #XMMEA_EM/P in 10 - 15 keV. As a result, we obtained good 
exposures, as shown in table 3. For the spectral analysis, MOS events with PATTERN < 12, 
#XMMEA_EM and FLAG = 0, and pn ones with PATTERN < 4, #XMMEA_EP and FLAG 
= 0, were finally selected. 

We integrated the EPIC spectra of NGC 2403 Src 3 and BGD, within a circle of 36" 
radius and an annulus of 72"-108" radius, respectively, both centered on Src 3. Because the 
pn image of the third observation revealed a faint point-like source within the BGD annulus, 
a circle of 30" radius around this contaminating source was removed from both MOS and pn 
data. We created the rmf and arf, utilizing the SAS tools, rmf gen and arf gen respectively. 

4-3. Spectral modeling 

Figure 4 compares the Chandra and XMM-Newton spectra of NGC 2403 Src 3 with the 
Suzaku one, in the form of the unfolded spectrum (panel a) and of the ratio to the Suzaku best- 
fit MCD model (panel b). All the data are presented in 0.7 - 7 keV, where the energy ranges of 
all the three instruments overlap. In the Chandra (4628) and XMM-Newton (0164560901) ob- 
servations, the spectral shape of Src 3 thus coincides well with those in the Suzaku observation, 
indicating that the source was in the MCD-like spectral state on these occasions. In addition, 
the spectral normalization is similar, to within ~20%, among these three data sets. In contrast, 
the X-ray spectrum in the Chandra observation of ObsID = 4630 exhibits a distinct and less 
convex shape, which is approximated by a power-law, at least above ~ 1 keV where the effect 
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of interstellar absorption become unimportant. The flux on this occasion is estimated to be 
lower by about 15% than that in the Suzaku observation. In figure 4, the deconvolved spectra 
have been derived using the absorbed MCD or PL model (whichever is more successfull) to be 
described below. 

The raw Chandra and XMM-Newton spectra of Src 3 in the individual observations 
are shown in figures 5 and 6, respectively. We applied the MCD and PL models to these 
spectra, in the same way as for the Suzaku data, and obtained the results summarised in table 

4. We neglected the contamination from the diffuse X-ray emission associated with the host 
galaxy, because the integration area, we adopted for Chandra and XMM-Newton data, were 
more than an order of magnitude smaller than that for Suzaku data. As figure 4 suggests, the 
MCD model better reproduced the observed X-ray spectra in general, except for that from the 
Chandra observation of ObsID = 4630. When the source was in the MCD-like state, the spectral 
parameters exhibited only small deviations from the Suzaku measurements. The bolometric 
luminosity of Src 3 was found to vary only about 10% among these observations. Similar MCD 
parameters were reported by Feng & Kaaret (2005) and Stobbart et al. (2006), based on the 
third XMM-Newton observation. 

The Chandra spectrum of ObsID = 4630, in contrast, was better reproduced by a PL 
model with a photon index of T = 2.37 ±0.08. The X-ray flux on this occasion, 9.3 x 10~ 13 
ergs cm -2 s _1 in the 0.7-7 keV range, was lowest among the data sets analyzed here. The 
absorption column density became larger by a factor of about 2-3 than those in the other 
observations; JVh = 3. 73^33 x 10 21 cm -2 . Similarly, an apparent change was reported from 
a transient ULX, CXOU J132518. 2-430304, which was discovered within the host galaxy of 
Centaurus A (Sivakoff et al. 2008). However, such significant changes in ./Vh, associated with 
transitions between the PL-like and MCD-like states, are uncommon among Galactic BHBs. 
Therefore, this Nn change from Src 3 is likely to be an artifact, caused presumably because 
a simple PL model is inappropriate to describe the softest end of the Chandra spectrum from 
ObsID=4630. We discuss this issue briefly in section 5. 

Although the MCD model gave a better fit than the PL model to all the XMM-Newton 
data, the third XMM-Newton dataset still exhibit some residuals with y 2 /d.o.i. = 637.9/568 
(figure 6c). This urged us to apply the variable-p disk model to this dataset and the other 
MCD-like ones, as we did for the Suzaku spectrum. Then, all the spectra in the MCD-like 
states were successfully reproduced by the variable-p disk model. The fit was significantly 
improved in the third XMM-Newton observation and the Chandra ObsID = 4628, and yielded 
p = 0.58±HI and 0.601HL respectively. 

5. Discussion 

In the Suzaku observation of NGC 2403 performed on 2006 March 16-17, Src 3 exhibited 
a curved 0.4 - 10 keV spectrum, which can be reproduced by the MCD model. We obtained 
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the innermost disk temperature and radius of the source as T in = 1.08_ a0 3 keV and R m = 
122.li 6 ; 8 a 1 / 2 km, respectively. The bolometric luminosity of the source was measured as 
Lboi = 1-82 x 10 39 a ergs s _1 . Making use of the archival Chandra and XMM-Newton data, 
we found the source to exhibit similar MCD-type spectrum on 4 occasions, while a PL-shaped 
spectrum on one occasion (Chandra observation of ObsID = 4630) when the 0.7-7 keV flux was 
about 15% and 17% lower than those in the Suzaku observation and the Chandra observation 
of ObsID = 4628, respectively. This flux change between the MCD-like and PL-like states 
are larger than the relative flux uncertainties between Suzaku, Chandra, and XMM-Newton 
(< 10%; e.g., Reeves et al. 2007; Bamba et al. 2008; Ishida et al. 2008). 

In the MCD-like states, T; n and Lboi of Src 3 exhibited only small variations across more 
than 8 years from the ASCA observation (Kotoku et al. 2000) to the Suzaku one. Both of 
these parameters are located at the lowest end of those of typical ULXs. The value of i?; n 
derived in these observations corresponds to the radius of the last stable orbit, 3Rs (Rg being 
the Schwarzschild radius), around a non-rotating BH with a mass of ~ 13a 1//2 M . Because 
the Eddington luminosity of such a BH becomes Le ~ 2 x 10 39 a 1//2 ergs s _1 , the properties 
of Src 3 in the MCD-like states can be understood consitently as a ~ 13M Q BH radiating at 
close to Le, wherein a standard accretion disk is formed down to ~ 3i?s, like in the high/soft 
state of normal Galactic BHBs. This apparently reconfirms the conclusion by Kotoku et al. 
(2000). Then, can we really explain the overall behavior of Src 3, by assuming that it was in 
the standard high/soft state during the Suzaku observation ? 

One difficulty with the high/soft state interpretation comes from the fact that the third 
XMM-Newton dataset preferred a variable-p disk model with p = 0.58^g;o2, which is smaller than 
that in the standard accretion disk. To see this in a systematic manner, we plotted in figure 7, 
the parameters determined from the variable-p disk model, against the observed X-ray flux in 
the 0.7-7 keV range. Although errors are rather large, the plot reveals a tendency that R in 
and p become smaller, and T m gets higher, as the flux increases. Such behavior is numerically 
predicted by the slim disk model (e.g., Watarai et al. 2000; Ohsuga 2006). The correlation is 
not affected significantly by possible flux uncertainties among the three sattelites, because the 
three XMM-Newton data points reveal the same behavior. Therefore, it is natural to regard 
that the source was actually in the slim disk state during these observations. Moerover, figure 7 
indicates that Src 3 was in the transition regime between the standard high/soft state (p = 0.75) 
and the slim disk regime (p < 0.75) when observed with Suzaku. 

Another difficulty with the high/soft interpretation is provided by the spectral transition, 
revealed by the Chandra data (ObsID = 4630), when the source was ~ 15% less luminous than 
in the Suzaku data. The spectral slope of this PL-like state, T = 2.37 ±0.08, is too steep 
for the source to be interpreted as in the low/hard state (r = 1.5 - 2.0; e.g, Tanaka 1997), 
which would be realized when the source becomes less luminous than in the high/soft state. In 
contrast, this steep slope agrees very well with those of Galactic BHBs in the VHS (Miyamoto 
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et al. 1991; Kubota & Done 2004). Furthermore, Galactic BHBs have been observed to become 
less luminous when they make a transition from the slim disk state to the VHS (Kubota & 
Makishima 2004; Abe et al. 2005). These comparisons suggest that the PL-like state of Src 3 
corresponds to the VHS. 

Since Comptonization of disk photons by hot coronae is thought to become important 
in the VHS (Kubota & Makishima 2004), we tried a Comptonization model, compTT in 
XSPEC (Titarchuk 1994), on the Chandra spectrum of Src 3 in ObsID = 4630. The electron 
temperature of the Comptonizing coronae was fixed at T e = 20 keV, as seen in the VHS of XTE 
J1550-564 (Kubota & Done 2004). The model successfully described the data with parameters 
shown in table 5, giving a fit goodness similar to that with the PL model. The optical depth of 
the corona, r = 1.2 ±0.1, is comparable to those of the Galactic BHB in the VHS (r ~ 2; e.g, 
Kubota et al. 2001; Kubota & Done 2004). Furthermore, the value of N-r has become consistent 
with those found in the MCD-like states; this is because the Comptonized continuum flattens 
in energies below that of the seed photons (T see< j ~ 0.25 keV). Although the derived value of 
T see d is lower than the disk temperature T m ~ 1.1 keV measured in the MCD-like state, the 
discrepancy could be due to a disk truncation at a radius larger than that of the last stable 
orbit, as is suggested in the "strong" VHS of XTE J1550-564 (Kubota & Done 2004). 

From these considerations, we conclude that NGC 2403 Src 3 normally resides at the 
lowest end of the slim disk state and makes occasional transitions into the VHS. Since the 
dominance of these two spectral states is regarded as an essential property of the ULXs (e.g., 
Kubota et al. 2002; Tsunoda et al. 2006; Mizuno et al. 2007), we may regarded Src 3 as a ULX, 
in spite of its moderate luminosity. In addition, we may regard Src 3 as radiating at ~ Le, for 
the following two reasons. One is that the transisions of BHBs between the VHS and the slim 
disk state is obvseved at > 0.3Lg (Kubota & Makishima 2004; Abe et al. 2005), while the other 
is that the value of p ~ 0.75 would not be observed if the source were radiating with ^> Lg, 
Thus, Lboi of the source gives the BH mass as M ~ 10 - 2OM . 

We may independently estimate the BH mass using i? in determined by the variable-p 
disk model. In this case, the simple assumption of i?; n = 3Rs would underestimate a mass of the 
BH that hosts a slim disk, because the X-ray photons would also be radiated from the region 
inside the last stable orbit (Watarai et al. 2000). Recently, Vierdayanti et al. (2008) evaluated 
a factor to correct the underestimation as 1.2 - 1.6, by analyzing numerically simulated X-ray 
spectra from the slim disk with the variable-p disk model (the extended disk black body model 
in their paper) for a wide range of mass accretion rate. Combining this correction factor and 
i?i n derived from our variable-p disk model fit (figure 7b), we estimate the BH mass of NGC 
2403 Src 3 as M = 9 - 15M Q . This value agrees very well with that derived from the Eddington 
limit argument. Therefore, we can consistently interpret the available data of NGC 2403 Src 
3, by considerig that it hosts a BH with a mass of M = 10 - 15M Q , radiating at ~ 

Finally, let us discuss implications of the present results on the ULX phenomenon in 
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general. Obviously, they do not give direct support to the intermediate mass BH interpretation 
of ULX, since the observed luminosity of Src 3 is well within the Eddington limits of ordinary- 
mass stellar BH. Nevertheless, one important discovery is that this relatively low-luminosity 
object exhibits clear ULX behavior, including in particular a transition between the slim disk 
state and the VHS. In other words, the ULX phenomenon appears even at a luminosity range 
of ~ 2 x 10 39 ergs s _1 . This is within ~ 3 times the Eddington luminosity of any stellar-mass 
BHs, even if we consider the least massive known BHBs as GRO J0422+32 (3.97 ± 0.95M Q ; 
Gelino & Harrison 2003), GRO J1655-40 (5.5-7.9M ; Shahbaz et al. 1999), and GX 339-4 
(5.8±0.5M Q ; Hynes et al. 2003). Consequently, the ULX behaviour in this case does not 
require any highly (~ 10) super-Eddington luminosity. Extrapolating this conclusion to more 
luminous ULXs, we suggest that they are not radiating at highly super-Eddington lumi- 
nosities, either, even though their mass accretion rates are very likely to be "super-critical". 
That is, ULXs are thought to be BHs which are significantly more massive than Galactic 
BHBs, accreting matter at super-critical rates, and yet radiating at about their Eddington limit. 

We are grateful to all the members of the Suzaku team, for the successful operation 
and calibration. We also thank Dr. Ohsuga for his valuable theoretical comments. This 
research has made use of the archival Chandra data and its related software provided by the 
Chandra X-ray Center (CXC). A part of the result is based on the observation obtained with 
XMM-Newton, an ESA science mission with instruments and contributions directly funded 
by ESA Member States and NASA. We have made extensive use of the NASA/IPAC Extra 
galactic Database (NED; the Jet Propulsion Laboratory, California Institute of Technology, the 
National Aeronautics and Space Administration). The optical DSS image of NGC 2403 was 
downloaded from SkyView 2 . The ROSAT source catalog were taken from the ROSAT Source 
Browser by Dr. Englhauser 3 . 
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Table 1. Best-fit parameters to the Suzaku spectra of Src 3. 



Model 


MCD 


variable-p 


PL 


N H (10 21 


cm- 2 ) 1.40t°;« 


1 co+l-02 
i - oz -0.74 


8.661°^ 


Ti„ (keV) 




1.09 ±0.07 


- 


p,T* 




0.73l°- 


2.79 ±0.07 


R m (km) 


i22.ii5;2 


116.4ll; 4 7 




L t 


1.82 


1.89 


4.10 


L th * 


0.035 


0.038 


0.174 


X 2 /d.o.f. 


182.0/174 


181.9/173 


349.4/174 



* Index p ot radial temperature dependence tor the variablc-p model, 
or the photon index T for the PL one. 



Lbol f° r the MCD and variable-p model, or absorption-corrected 0.5 
- 10 keV luminosity for the PL model, both in 10 39 ergs s . 
•f Absorption-corrected 0.5 - 10 keV luminosity of the sum of the two 
APEC components, in 10 39 ergs s^ 1 . 

Table 2. Log of archival Chandra observations of NGC 2403, utilized in the present paper. 



ObsID 


4628 


4630 


Date 


2004 Aug 


23 2004 Dec 22 


Exposure (ks) 


46.5 


48.7 


Mode 


Timed Exposure 


Frame Time 




3.2 s 


Format 




FAINT 


ACIS Chip* 


S3 


S2 


A6>t 


5'.8 


4'. 9 




15".7 


13". 8 


^BGD § 


19". 6-23" 


.6 17". 7-21". 6 


Signal (cts s" 1 ) 11 


0.16 


0.11 



* ACIS CCU chip on which Src 3 is located, 
t Off axis angle of Src 3 from the aim point 

* Radius of the source integration circle. 

§ Radius of the BCD integration annulus. 
I Count rate of Src 3 in 0.3 - 10 keV 
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Table 3. Log of archival XMM-Newton observations of NGC 2403. 



\J DS1U 




ni ^n^^i i ni 


U10UD01ZU1 


ni f\A ^noni 


uate 




zUUo Apr oU 


zUUo oep ll 


zuu4 aep iz— io 


IVlUCic 






run rifiiiic 




Filter 




Thin 


Thin 


Medium 


Exposure (ks) 


MOS1 


5.26 


6.98 


57.1 




MOS2 


5.52 


7.20 


56.2 




pn 






49.6 



Table 4. Summary of fitting to the Chandra and XMM-Newton spectra of Src 3 



Instrument ObsID 


Model 


7V H 

(10 21 cm" 2 ) 


Tin 

(keV) 


p,T* 


-Rin 

(km) 


L t 


X 2 /d-o.f. 


Chandra 4628 


MCD 


1 47+0.15 

-L- 4 ' -0.14 


1.06 ±0.04 




127.4±|;| 


1.86 


166.9/156 




variablc-p 


2.23 ±0.44 


1 Oq + 0.15 

1 - ZO -0.11 


u - DU -0.05 


70.3l 2 «i 


2.48 


158.4/155 




PL 


q oq+0.26 
°- OO -0.24 




9 90+O.O7 
z - zo -0.06 




2.62 


230.0/156 


4630 


MCD 


Q7+ - 21 
u - y ' -0.19 


1 04+0-05 




115-3™ 


1.42 


194.8/123 




variable-p 


q rc+0.34 
J - oo -0.32 


3.3(>2.3) 


0.46 ±0.01 


4.193$ 


10.4 


109.9/122 




PL 


q 7O+0.35 

J - ' °-0.33 




2.37 ±0.08 




2.14 


110.7/123 


XMM-Newton 0150651101 


MCD 


1 40+ - 43 

1 -* u -0.38 


1.09 ±0.10 






1.72 


31.1/37 




variable-p 


1 44+1-31 


1 io +0 - 34 


n 74+I 8I 

u - '^-0.18 


ii2-4^; 9 8 


1.75 


31.1/36 




PL 


q 77+0.76 
°- ' ' -0.66 




9 9 q + 0.18 

z - zo -0.17 




2.42 


46.3/37 


0150651201 


MCD 


1 66 +0 ' 41 


1.11 ±0.09 




n6-5lf 7 ; 8 6 


1.87 


33.7/45 




variablc-p 


9 14+1-15 
^• 1 ^-1.22 


1 91+0.38 

1 - zl -0.23 


0-651°;?° 


oq 1 +110.4 
S ' ;i - i -48.5 


2.16 


33.2/44 




PL 


4 - lo -0.64 




9 99+O.I6 
z - zz -0.15 




2.65 


47.5/45 


0164560901 


MCD 


1.66 ±0.07 


1.06 ±0.02 




130.9ite 


1.95 


637.9/568 




variable-p 


2 5Q+ - 22 


1.27 ±0.06 


u -°°-0.02 


64 

°*-°-10.1 


2.85 


585.6/567 




PL 


4 S0+ 013 




2.43 ±0.03 




2.81 


943.3/568 



* Index p of radial temperature dependence for the variable-p model, or the photon index V for the PL model. 

T Lbol for the MCD and variable-p model, or absorption-corrected 0.5 - 10 keV luminosity for the PL model, both in 10 39 ergs s — 1 . 
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Table 5. Parameters of compTT fitting to the Chandra spectrum of ObsID = 4630. 



Parameters Value 



iV H (10 21 cm- 2 ) 


1 4+ 2 - 6 


T sood * (kcV) 


0.25(<0.29) 


T c t(keV) 


20 (fix) 


r* 


1.15 ± 0.11 



X 2 /d.oi. 111.14/122 

* Energy oi seed soft photons 

t Electron temperature of the Comptonizing corona 

* Optical depth of the Comptonizing corona 
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Fig. 1. Suzaku XIS contour image of NGC 2403 in the 0.5 - 10 keV range, overlaid on the optical DSS 
gray-scale image. Data from all the XIS CCD chip are summed up. Neither background-subtraction nor 
exposure-correction is performed to the image. The source signals of Src 3 and the associated background 
are integrated within Src3 and BGD3, while those of Src 5 are extracted from Src5 and BGD5. The 
positions of the ROSAT X-ray sources, taken from the second ROSAT source catalog of pointed obser- 
vations with the position sensitive proportional counter (2RXP), are indicated with crosses, and that of 
SN2004dj with the filled star. 
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BGD-subtracted Lightcurve of NGC 2403 Src3 

X|S Fl Bin time: 5760. s 
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-O- 
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2x10 4 4x10 4 6x10 4 8x10 4 
Time (s) 

Start Time 1 381 1 5:07:1 5:1 84 Stop Time 1 381 1 1 3:31 :1 5:1 84 

Fig. 2. Background-subtracted XIS FI lightcurve of NGC 2403 Src 3. Each time bin is set to 5760 s, 
corresponding to the orbital period of Suzaku. Panels (a) and (b) show the soft (0.5 - 2 keV) and hard (2 
10 keV) energy band lightcurves, respectively. The hardness, simply calculated as the ratio of the hard 
band count rate to the soft band one, is shown in panel (c). Dashed and dotted lines in panel (c) indicate 
predictions by MCD models with of T in = 1.0 keV and 1.2 keV, respectively. 
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NGC 2403 Src 3 

Suzaku Observation 
0.1 „ , ,- 




Fig. 3. (a) Suzaku XIS spectra of NGC 2403 Src 3, shown together with the best-fit MCD model. The FI 
and BI data are shown with black and red points, respectively. Two APEC components are also indicated 
(see text). The data are binned into pixels with at least 100 events. Panels (b) and (c) show the residuals 
from the MCD and PL models, respectively. 
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NGC 2403 Source 3 
(a) Unfolded Spectra 
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Fig. 4. (a) Unfolded spectra of NGC 2403 Src 3. The black, red, green and blue points indicate the 
Suzaku data, the Chandra data of ObsID = 4628, those of ObsID = 4630, and the XMM-Newton data of 
ObsID = 0164560901, respectively. Panels (b) - (d) show the Chandra and XMM-Newton spectra, divided 
by the prediction of the best-fit MCD model determined by the Suzaku spectrum. The color specification 
is the same as in panel (a) 
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NGC 2403 Source 3 NGC 2403 Source 3 

(a) ObsID = 4628 (b) ObsID = 4630 




0.5 1 2 5 0.5 1 2 5 

Energy [keV] Energy [keV] 



Fig. 5. Chandra ACIS spectra of NGC 2403 Src 3, from ObsID = 4628 (panel a) and ObsID = 4630 
(panel b). The employed model is indicated in each panel. 
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NGC 2403 Source 3 
(a) ObslD = 01 50651 101 
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NGC 2403 Source 3 
(b) ObsID = 0150651201 
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Fig. 6. XMM-Newton EPIC spectra of NGC 2403 §£g 3, obtained on 3 occasions. The black, and red data 
points indicate the MOS and pn spectra, respectively. The histogram in the individual spectra indicate the 
best-fit MCD model. For the observation of ObsID = 0164560901 (panel c), the residuals of the variable-p 
disk model are also plotted. 
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Fig. 7. A summary plot of the variable-p modeling of the Suzaku, Chandra, and XMM-Newton spectra 
of Src 3. Panels (a), (b) and (c) present p, R[ n and T| n as a function of the observed 0.7 - 7 keV flux, 
respectively. The dotted line in panel (a) indicates p = 0.75. 
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A. 0.1. X-ray spectra of NGC 2403 Src 5 

The X-ray signals from NGC 2403 Src 5 is clearly detected in figure 1. We breifly analyze 
the X-ray spectrum of Src 5, altough we regard that a detailed study of the source is beyond 
the scope of the present paper, due to relativley low signal statistics. 

The X-ray signals from Src 5 and background are extracted from the circles denoted 
as Src5 and BGD5 in figure 1, respectively. Figure 8 shows the XIS spectra of Src 5. After 
calculating the rmf and arf files in the same way as for Src 3, the MCD and PL models were 
applied to the spectrum. We did not include the thermal components from the galaxy in the 
fitting, since the signal statistics are rather low and Src 5 is located slightly outside the diffuse 
emission (see figure 1). The absorption column density for the MCD model was fixed at the 
Galactic value, since it was unconstrained from the data. 

The results of this analysis are summarised in table 6. Thus, the MCD model was mod- 
erately successful, and yielded the disk parameters as T m = 1.52taii keV and i?i n = 28. 6i 3 ; 7 a 1//2 
km. The bolometric luminosity was estimated to be L^ Q \ =4.0 x 10 38 a ergs s" 1 . The PL fit was 
slightly less successful, and gave a photon index of Y = 1.65lo]i2- 

Table 6. Best-fit parameters for to Suzaku spectra of Src 5. 





MCD 


PL 


iV H (10 21 cm- 2 ) 


0.41+ 


1 1O+0-73 


Tin (kcV) 


1 r 9 + 0.13 




R in (km) 


28.6±f;£ 




r 




1 fi5+ 013 


L * 


4.0 


4.7 


X 2 /d.oi. 


138.8/121 


150.5/121 



* ^bol t° r the M(JL) model, or absorption- 
corrected 0.5 - 10 keV luminosity for the PL 
ones, in 10 38 ergs s _1 . 
i fixed at the Galactic value. 
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NGC 2403 Src 5 
Suzaku Observation 
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